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Abstract. We present visible, near infrared, and mm-wave observations of RCW 108, a molecular cloud complex 
in the Ara OBI association that is being eroded by the energetic radiation of two 0-type stars in the nearby 
cluster NGC 6193. The western part of the RCW108 molecular cloud, for which we derive a mass of 8000 M0, 
contains an embedded compact HII region, IRAS 16362-4845, ionized by an aggregate of early-type stars for which 
we estimate a mass of ~ 210 Mq. The spectral type of the earliest star is 09, as confirmed by the visible spectrum 
of the compact HII region. We notice a lack of stars later than AO in the aggregate, at least having the moderate 
reddenings that are common among its B-type stars, and we speculate that this might be a consequence of the 
extreme youth of the aggregate. We also note the existence of a dense ionized clump (n > 10^ cm~^) appearing 
near the main ionizing star of the compact HII region. We examine the distribution of stars displaying infrared 
excesses projected across the molecular cloud. While many of them are located in the densest (n ~ 10*"^ cm~^) 
area of the molecular cloud near the position of IRAS 16362-4845, we also find a group concentrating towards 
the edge of the cloud that faces NGC 6193, as well as some other stars beyond the edge of the molecular cloud. 
The intense ionizing radiation field by the O stars in NGC 6193 is a clear candidate trigger of star formation in 
the molecular cloud, and we suggest that the existence and arrangement of stars in this region of the molecular 
cloud supports a scenario in which their formation may be a consequence of this. However, infrared excess stars 
are also present in some areas of the opposite side of the cloud, where no obvious candidate external trigger is 
identified. The existence of such tracers of recent star formation scattered across the more massive molecular 
cloud associated with IRAS 16362-4845, and the low star formation efficiency that we derive, indicate that it is 
in a state to still form stars. This is in contrast to the less massive cloud (~ 660 M0) close to NGC 6193, which 
seems to be more evolved and mostly already recycled into stars, and whose internal kinematics show hints of 
having been perturbed by the presence of the massive stars formed out of it. 
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1. Introduction 

RCW 108 (Rodgers et al. is an extended HII re- 

gion belonging to the Ara OBI association (Herbst & 
Havlen nWT!, Kaltcheva & Georgiev '1992?, highhghting 
one of the few remaining places of Ara OBI where star 
formation is still active (Yamaguchi et al. I1999|l . The ex- 
Send offprint requests to: F. Comeron 

* Based on observations obtained at the European Southern 
Observatory using the ESO New Technology Telescope (NTT) 
(programs 61.D-0447, 64.L-0049, and 71.C-0429), the ESO- 
MPI 2.2m telescope (program 62.1-0454), the ESO 3.6m 
telescope (program 69.C-0522(A)) and the Swedish-ESO 
Submillimeter Telescope (program 61.C-0243), on La Silla, 
Chile. 



istence and arrangement of objects and structures usu- 
ally associated to massive star formation (early-type stars, 
molecular clouds with dense cores, embedded pre-main 
sequence stars, and ionized gas) makes it an interesting 
target for the study of the interplay between the inter- 
stellar gas and dust and newly formed stars. Its distance 
of approximately 1.3 kpc (see Arnal et al. (2003) for a 
brief discussion; this is also the value that we adopt in 
the present paper) allows a study at good spatial reso- 
lution (1' equals -0.4 pc). The HII region, NGC 6188, 
is actually the bright rim of a molecular cloud containing 
several star formation sites. The eastern edge of this cloud 
is being eroded by the ionizing radiation of the nearby O- 
type stars HD 150135 and HD 150136 in the open cluster 
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NGC 6193. Embedded in this cloud lies a compact HII 
region (Shaver & Goss 1197011 coincident with the source 
IRAS 16362-4845. This region is easily noticed in visible 
images of the area, but it is most prominent at near in- 
frared wavelengths at which the diffuse bright-rimmed HII 
region is hardly seen. The study of the ionizing stars of the 
compact HII region is especially interesting in the context 
of recent ideas on the determinant role that Trapezium- 
like clusters may play on the birth of high mass stars (e.g. 
BonneU et al.EMIl BonneU & Bate 2002). 

A pioneering study of IRAS 16362-4845 was published 
by Straw et al. (19871), who presented infrared imaging 
observations ranging from the J band (1.25 ^m) up to 
100 fjLVH, as well as low resolution infrared spectroscopy of 
two selected sources and of the ionized nebula. Although 
that work established some of the main characteristics of 
the compact HII region and its associated stellar popu- 
lation, its depth and spatial resolution are rather modest 
by the standards of current instrumentation on medium- 
sized telescopes. Moreover, only IRAS 16362-4845 and its 
immediate vicinity are considered in Straw et al. H1987|l . 
A recent study based on radio recombination line, ra- 
dio continuum, molecular line, mid-infrared (MSX), and 
near-infrared (2MASS) observations of a region compara- 
ble to the one discussed in the present paper has been 
published by Urquhart et al. I|2()()4|l . Finally, low angu- 
lar resolution molecular observations of RCW 108 in the 
13C0 J=1^0 line (Yamaguchi et al. lTM5|l and the ^^CO 
J=1^0 line (Arnal et al. 2003 ) provide useful complemen- 
tary information on the larger scale distribution of lower 
density gas in the region and the existence of other star 
forming sites in Ara OBI. 

This paper presents new visible and near infrared 
imaging and low resolution spectroscopy of IRAS 16362- 
4845 and its surroundings, aiming at complementing and 
updating the analysis of Straw et al. H1987|l and at ex- 
tending it to a broader area of the molecular cloud. We 
also present maps in the lines of ^^CO J=2— ^1 and ^■^CO 
J=l— >-0 centered on IRAS 16362-4845 covering most of 
the area included in our near infrared images, and south of 
the emerged cluster NGC 6193. Additional interferometric 
observations of the Ha emission throughout much of the 
region centered on IRAS 16362-4845 are shown and dis- 
cussed. Our imaging and molecular-line observations thus 
cover the densest areas of the molecular clouds, whereas 
the low resolution spectroscopy focuses on the compact 
HII region and its embedded stellar component. 

2. Observations 

2.1. Near infrared imaging 

A near infrared mosaic in the J (1.25 /im), H (1.65 /im), 
and Ks (2.2 /im) bands, covering an area of 13' x 13' on 
the sky, was obtained on the night of 22/23 February 2000 
using the Sofl infrared spectrograph and array camera at 
the ESO New Technology (NTT) telescope. The central 
5' X 5' area containing IRAS 16362-4845 was imaged at a 



greater depth by means of stack of 15 frames in each filter, 
each one containing 6 individual exposures of 2 sec coad- 
ded on the detector, obtained with small telescope offsets 
in between. The peripheral area was imaged by pointing 
the telescope at 32 regularly spaced positions defining the 
sides of a square centered on IRAS 16362-4845, each frame 
being in turn the result of 6 individual exposures of 2 sec- 
onds coadded on the detector. The spacing between con- 
secutive pointings in the peripheral area was 1', i.e., 1/5 of 
the field of view of the camera. Each sky position was thus 
imaged by five different pointings, except near the borders 
of the square pattern. After flat fielding and dark subtrac- 
tion, the sky frame was constructed by median averaging 
the stack of the periphery frames (which were found to be 
virtually devoid of nebulosity) with clipping of the upper 
half of the pixel values. This sky frame was then sub- 
tracted from each individual pointing in both the periph- 
ery and the central area. All the ifg-band images were 
then registered to construct the mosaic, using the posi- 
tions of stars common to consecutive images as references 
to compensate for the telescope offsets. The mosaics in 
the other two filters were constructed by individually reg- 
istering each of their component frames with the i^s-band 
mosaic, in order to avoid slight relative scale distortions 
between filters due to accumulation of small errors in the 
registering process. 

Sources were detected in our frames using DAOFIND 
fStetson ll987|l . Relatively isolated unsaturated images of 
bright stars were used to determine an approximate PSF 
needed for the identification of point sources. Due to the 
crowdedness of the field outside the areas of densest nebu- 
losity and to our interest in point sources only, photometry 
on the resulting mosaic was performed by defining an un- 
dersized aperture at the position of each detected star, 
measuring the flux inside it, and then adding the rest of 
the flux in the PSF as given by the fit of a circularly sym- 
metric radial profile to each stellar image. This procedure 
allowed both to remove the contamination to the photom- 
etry to other stars located on the wings of the PSF, and 
to adjust to the mildly variable image quality across the 
field of view. 

2.2. Visible spectroscopy 

Spectroscopy was carried out on the night of 1/2 April 
2003 using EMMl, the visible imager and spectrograph 
at the NTT. We chose a grism yielding a coverage of 
the 3800 A < A < 9700 A interval at a resolution of 
A/ A A = 570 with the l"0-wide sht that we used. The slit 
was placed on the line joining the brightest star seen pro- 
jected on IRAS 16362-4845 in visible-light images and the 
faint, very red star 4"0 to its Southwest that becomes the 
brightest one in the near infrared K band^. The effective 
slit length, 8'0, provided us with a cut across the com- 

^ These are respectively denominated Star 8 and Star 12 
in our discussion of the stellar aggregate associated to 
IRAS 16362-4845 in Section 1^ 
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pact HII region but also included a bright segment of the 
rim nebula. We obtained 3 separate spectra of 10 minutes 
of exposure time each on the same telescope position, and 
then coadded the resulting frames. Spectra of the two stars 
were extracted using the APALL task on IRAF, defining 
sections adjacent to the traces of the stars for the subtrac- 
tion of the sky and nebular emission together. Spectra 
of the nebula at selected points were extracted as well 
by interactively defining the apertures on the frames and 
defining an appropriate sky aperture in a suitable area 
(see discussion in Section f3.5|l . Relative flux calibration 
was performed using the spectrum of LTT 4364 (Hamuy 
et al. 11992(1 as a reference. Wavelength calibration was 
carried out by extracting the spectrum of a ThAr lamp 
obtained with the same instrumental setup at the same 
aperture positions. 

2.3. Near infrared spectroscopy 

Near infrared simultaneous low resolution spectroscopy 
in the H and K bands was obtained using SofI in spec- 
troscopy mode on the night of 2/3 April 2003. The grating 
used covered the 1.5 fim - 2.4 yum interval at a resolu- 
tion A/ A A — 590 with our 1"0 slit. In this case the slit 
was placed on the line joining the brightest near infrared 
source and another source 3" 5 to its Southwest that is sim- 
ilarly bright in the infrared but redder in color, thus being 
undetected in our visible images^. Four individual spectra 
of 50 sec exposure each were obtained by slightly offsetting 
the telescope between exposures along the direction of the 
slit. Spectra were then extracted using a procedure anal- 
ogous to that followed for the visible spectra. Cancelation 
of telluric features in the extracted spectra was achieved 
by observing the G3V star HIP 81746 immediately after 
the observation of RCW 108 at a very similar airmass, and 
relative flux calibration was carried out by assuming that 
the overall shape of the spectrum of HIP 81746 is well ap- 
proximated by a black body at a temperature of 5700 K 
over the 1.5 fim - 2.4 /im interval. Wavelength calibration 
in the infrared was performed by using the airglow OH 
emission lines as a reference (Oliva & Origlia|TO92,). 

2.4. Additional imaging observations 

Besides the imaging and spectroscopy described in the pre- 
vious Sections, we obtained some additional visible imag- 
ing that is helpful in illustrating the overall, large scale 
morphology of RCW 108 as well as the more detailed 
structure of the IRAS 16362-4845 HII region. For this pur- 
pose, a Ha image of the area was obtained using the Wide 
Field Imager (WFI) at the ESO-MPI 2.2m telescope on 
La Silla on the night of 26/27 March 1999. The image 
was obtained by combining four individual pointings of 
the telescope with offsets of approximately 2' in order to 

^ These are respectively denominated Star 12 and Star 16 
in our discussion of the stellar aggregate associated to 
IRAS 16362-4845 in Section 1^ 



cover the gaps between the individual CCD chips com- 
posing the WFI detector array. The exposure time of each 
individual frame was 300 sec. Additional short exposures 
through the B and V filters were also obtained, allowing 
their combination into a color image that has been widely 
reproduced elsewhere (see e.g. Collins Petersen I2001f) . 

We also obtained U BVRI images of a smaller, 5'3 x 5'3 
field centered on IRAS 16362-4845 using the SUSI2 visible 
imager on the NTT on the night of 26/27 August 1998. 
Since the SUSI2 detector is also composed of two chips 
separated by a gap we divided the exposures into three 
pointings separated by small telescope offsets, allowing 
the reconstruction of a continuous image of the field with 
the gaps filled. The total exposure time in each filter is 9 
minutes. 

2.5. IVIillimeter CO observations 

Observations of ^^CO J=2^1 at 230.538 GHz and ^^CO 
J=1^0 at 110.201 GHz were carried out using the 
SEST during September 1998 towards IRAS 16362-4845 
(an area of ~ 10' x 8' was mapped) and south of the 
cluster NGC 6193 (4' x 7'). We refer hereafter to the 
molecular gas detected at those positions as the west- 
ern and the eastern clouds, respectively. The map cen- 
ters were a(2000)=16''40"00n3, (5(2000)=-48°51'46.5" 
and a(2000)=16''41"20"32, (5(2000)=-48°47'45.1" re- 
spectively. A pointing grid of 22 arcsec was adopted, sim- 
ilar to the beamsize at 230 GHz (23") and half-beamsize 
at 110 GHz (45"). The velocity resolution of all data is 
0.11 kms~^ and the velocity coverage is —32 to —7 km s~^ 
for the eastern cloud (around NGC 6193) and —36 to 
— 13 km s^^ for the western cloud (around IRAS 16362- 
4845). The average system temperature throughout the 
observations was 210 K for ^^CO J=2^1 and 140 K for 
^^CO J=1^0 . The calibration was done using the stan- 
dard chopper-wheel method. The resulting atmosphere- 
corrected antenna temperatures were converted to main- 
beam brightness temperatures using the values for main- 
beam efficiencies quoted in the SEST Handbook (0.5 for 
230 GHz and 0.70 for 110 GHz). The telescope pointing 
and subreflector focusing were checked regularly, we esti- 
mate the pointing accuracy to be better than 5 arcsec and 
adopt the standard SEST value of 10% for the uncertainty 
in the antenna temperature scale. The average r.m.s. main 
beam brightness temperature noise per channel is 0.5 K 
for i^CO J=2^1 and 0.2 K for ^^CO J=1^0 . 

2.6. Interferometric Ha observations 

The Ha interferometric observations were made with the 
CIGALE instrument on the 3.6m telescope (La Silla) in 
April 2002. The data cubes were obtained with a spa- 
tial resolution of 0.4" for 4 fields centered on the com- 
pact HII region. The size of each field is 4.5 arcmin. The 
Fabry-Perot interferometer used has an interference order 
1938 (at Ha wavelength) providing a spectral sampling 
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Fig. 1. A negative Ha wide-field image of the RCW 108 area fSect. inji . showing the overall distribution of the ionized 
gas in the region and main objects discussed in the text. The large-scale distribution of the column density of obscuring 
dust can be inferred as well by noting the highly variable surface density of background stars across the field. The 
rings appearing around or next to bright stars are out-of-focus ghost images caused by internal reflections. 



of 3.2 km and a free spectral range of 155 km s~^. 
The interference filter used is centered at 6562 A with a 
FWHM of 11 A. The velocity and line width accuracy is 
estimated to be 0.7 km s~^. A complete description of the 
instrument, including data acquisition and reduction tech- 
niques has been given in le Coarer et al. |I992). The Ha 
profiles are decomposed with the Ha geocoronal night-sky 
line (modeled by a purely instrumental profile) and a neb- 
ular line (modeled by an instrumental profile convolved 
with a gaussian). An automatic procedure has been used 
to perform the fit. Each field contains 512x512 pixels, but 



in order to increase the Signal-to-Noise ratio (S/N) we ex- 
tracted and analysed profiles from areas of 6.5" x 6.5" size. 
A rough flux calibration was performed by observing the 
planetary nebulae NGC 2899 and adopting fluxes given 
by Perinotto and Corradi 1*11998,1 . 
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3. Results 

3.1. Morphology of the distribution of gas 

The overall distribution of ionized and molecular gas over 
the RCW 108 area and its surroundings is well illustrated 
by the negative wide-field Ha image presented in Figure^ 
The cluster NGC 6193 occupies the central/eastern side 
of the field and is dominated by the close pair of O stars 
HD 150135/150136. Glow in Ha - visible as dark patches - 
is present over most of the field, mostly arising from gas on 
the background of NGC 6193, as well as probably from gas 
left over from the formation of the cluster indicated by a 
peak in molecular line emission seen in the same direction 
as its brightest stars (see Section lT^ . A simple inspection 
of the starcount density in the area to the East of the rim 
nebula indicates that the extinction on the background is 
low in the Northern half and higher towards the Southeast. 

The arrangement of illuminated and shadowed areas 
in the rim nebula suggests that most of the interface be- 
tween the molecular cloud and the ionization front eroding 
it is seen roughly edge-on. The pattern of bright nebulosi- 
ties and shadows dramatically reveals the intricate three- 
dimensional structure, as well as the existence of regions 
with widely varying densities. The column density of the 
gas near the Northern and Southern edges of the rim neb- 
ula is low, as seen from the only slight decrease of stellar 
density in its direction, as compared to the more opaque 
clouds located near the center of the image. Small, high 
density cores are visible as dark patches projected against 
a brighter background all over the nebula, and some of 
them appear on areas where the erosion front has trav- 
eled past them, leaving them isolated from the bulk of the 
molecular gas to the West. 

The most conspicuous star forming site of RCW 108, 
IRAS 16362-4845, roughly coincides with the thickest part 
of the molecular cloud, as shown by both starcounts and 
molecular-line maps fSection l3.2l see also Figure|2J. North 
of it one finds the reflection nebula surrounding the B3V 
star CD-48 11039 (Herbst ITTFT^ . No other star forming 
sites are obvious in either visible or infrared images, al- 
though we will discuss evidence for other lower-mass star 
forming sites spread across RCW 108 in Section 

A comparison of visible and near infrared images of 
IRAS 16362-4845 and its surroimdings, shown in Figure^ 
provides useful information on the structure of the com- 
pact HII region. In the visible, the brightest patch of neb- 
ulosity lies at the position of three bright and only lightly 
reddened stars evenly spaced roughly in the North-South 
direction. The peak of visible Ha emission is displaced to 
the South of the nominal IRAS position, where the visible 
Ha emission is comparatively faint. This indicates that 
the core of the compact HII region is heavily reddened, 
and thus that the northernmost of the three bright visible 
stars is not its ionizing star, as already noted by Straw et 
aL itmTll . 

Broad-band visible images in B and V show reflection 
nebulosity to the Southeast of the zone of most intense 



Ha emission. This might be caused by a relatively unim- 
peded line of sight from that vantage point to the stars 
ionizing the HII region. However, we believe it more likely 
that the source of illumination is actually an anonymous 
star located to the Southeast, at a(2000) = 16''40'"08"21, 
J(2000) = -48°53'49"9. Although this star does not show 
any distinctive signs of belonging to the association from 
the data at hand, indirect evidence comes from a very 
red source lying 28" to its Northeast, which seems to be 
still embedded in the remnants of the core from which it 
formed. This core casts a shadow on the surrounding neb- 
ula that points directly away from the star noted above. 

The relatively unobscured view into the nebula pro- 
vided by the infrared images shows important differences 
with the visible-light picture. The K-hand image, which 
mainly traces the emission in Br7 (2.166 /im) and Hel 
(2.058 /im), is now clearly peaked at the position of the 
IRAS source, where a tight cluster of reddened stars 
(Section I3.4|l is seen, and extends eastwards in the gen- 
eral direction of NGC 6193. This eastward extension is 
totally blocked from view in visible images by an opaque 
layer of dust that in some small areas is thick enough to 
block the background emission even in the K band. 

To summarize, the overall structure of IRAS 16362- 
4845 hinted by the visible and infrared images is that of 
an embedded compact HII region. The foreground extinc- 
tion decreases towards the South, where the peak of the 
visible emission is, but the HII region itself presents an ex- 
tension towards the East that lies behind a thicker layer 
of obscuring dust. At the adopted distance of 1300 pc the 
size of the compact core is 0.16 pc (25"), surrounded by a 
fainter halo of 0.22 pc (35") in diameter, and the eastward 
extension reaches up to ~ 0.65 pc (100") from the center 
of the core. In the following section, we present the results 
of our molecular line mapping which support the scenario 
described above. 

3.2. Molecular line maps 
3.2.1. Overview 

We focussed on mapping the environment of IRAS 16362- 
4845 and NGC 6193 in the ^"^CO J==2^1 and 
"CO J=1^0 lines. The normally optically thin "CO 
J=l— >0 line is used to determine the column density and 
mass of the molecular clouds (Sec. 13. 2.311 while the ^^CO 
J=2^1 line provides an overview of the general small 
scale (23") distribution of molecular gas even on a low in- 
tensity level. Figure|21 (center) presents an overlay of ^^CO 
J=2^1 emission on the optical image where the eastern 
and western molecular clouds apparently associated with 
IRAS 16362-4845 and NGC 6193 are outhned by contours 
of CO emission. The western cloud corresponds very well 
to a region of high extinction (white areas in the Ha im- 
age) which is not the case for the eastern cloud fragment. 
There, regions of high extinction arc more diffuse and less 
clearly defined. 
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Fig. 2. A sequence of blue to infrared images of IRAS16362-4845 field, showing the changing appearance of the 
nebulosity as the extinction decreases when going to longer wavelengths. Visible images show all-pervasive, mostly 
reflection nebulosity, at least part of which is illuminated by stars not related to the HII region. The latter gains 
prominence in the red and infrared, revealing the aggregate of stars responsible for its ionization. 



Arnal et al. (200'6) mapped a 2.75° x3° region around 
RCW 108 at 8.7 resolution in ^^CO 1^0 and revealed 
a network of molecular clouds at velocities between —50 
and +6 km s^^. However, by using the galactic rotation 
curve they conclude that only CO emission between —27 
and —15 km s"^ is related to the Ara OBI association. 
This velocity range matches to our maps and indeed, we 
observe a very close correspondence between the western 
cloud and features related to Ara OBI (the rim nebula 
NGC 6188 and IRAS 16362-4845): the interface between 
HII region and western molecular cloud is impressively 
outlined by a sharp gradient of molecular line emission 
and follows closely the optical features. This indicates also 
that we see this region edge-on as already concluded from 
the Ha image alone (^Section l3.1|l . Another prominent fea- 
ture of the western cloud is the marked peak of ^^CO 
J=2^1 emission close to the position of IRAS 16362- 
4845. The ionized gas of this compact HII region is seen 
as a dark patch slightly shifted south-east from the main 
CO peak. Considering the distortion of the CO contour 
lines in this area, a close interaction between ionized and 



molecular gas is likely, i.e. erosion of molecular clumps by 
streaming ionized gas. 

In contrast to this rather evident correspondence be- 
tween ionized and molecular gas, the molecular emission 
close to the pair of O stars HD 150135/150136 does not 
show a clear morphology which is related to these stars. 
It shows a globular structure with a region of high col- 
umn density pointing away from the stars in south-east 
direction. Lower density material without directed struc- 
ture is present in the south and southwest. This cloud 
was mapped at 2' angular resolution by Phillips et al. 
H1986|l and they concluded that it represents a wind- 
swept globule probably shaped by NGC 6193, i.e. HD 
150135/150136. In our velocity integrated map, this sce- 
nario is not so clear and we will come back to this point 
in the next section where channel maps are discussed. 

A first impression of the complex velocity structure of 
the molecular clouds is given by the ^^CO J=2^1 and 
^•^CO J=1^0 spectra displayed in Fig. 13 The upper pan- 
els from the NGC 6193 region show that there are several 
velocity components which are present at all positions but 
vary in intensity. The most distinct feature is a single line 
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Fig. 3. The same negative Ha wide-field image of RCW 108 as shown in Fig. ^ is overlaid with contours of velocity 
integrated ^^CO J=2^1 emission at 23" angular resolution. The velocity range for the eastern (western) cloud is —24 
to -12 km s~i (-30 to -18 km s^^) and contours go from 17 (25) to 119 (300) K km s'^ in steps of 17 (25) K km s'^ 
Representative ^^CO J=2^1 (black line) and ^'^CO J=1^0 (filled grey) spectra taken at different positions from the 
eastern and western molecular clouds are shown in the panels above and below. For each cloud a positionally averaged 
^^CO J=2— >1 spectrum (grey line), amplified by a factor 4, is displayed as well. 



at —23 km s^^ whereas the other lines between —21 and 
— 12 km s^^ partly overlap. However, there is a clear veloc- 
ity gradient across the cloud, visible as a line shift com- 
pared with the the average ^^CO J=2— s-l spectrum (in 
grey) from —19 km (a) to —15 km s^^ (c). Since ^^CO 



J=2^1 and ^'^CO J=1^0 show the same line profile, ef- 
fects like self-absorption can be excluded in first order. 

The lower panels represent three positions from the 
western cloud where the line profiles are less variable. The 
^^CO J=2— >1 spectrum from the interface region (d) dis- 
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Fig. 4. -'^^CO J=2-^l (top) and ^''CO J=1^0 (bottom) emission distribution aroimd NGC 6193 in the velocity range 
-24.1 to -12.2 km s'^ (-20.3 to -13.5 km s'^). The velocity coverage of each ^^CO J=2-^l ("CO J=1^0 ) map 
is ~0.5 (0.3) km (five (three) velocity channels are comprised) and the central velocity is indicated in each panel. 
The stars represent HD 150135 and 150136. Contours go from 3a (1.5 K km s~^) to 453cr in steps of 30cr for -^^CO 
J=2^1 and 3a (0.4 K km s'^) to lOla in steps of 15cr for "CO J=l->0 . 



plays a flat-top profile (while the ^''CO J=1^0 line looks 
like a blending of several gaussian components) , probably 
indicating self-absorption effects at all velocities. In con- 
trast, the ^^CO J=2— >1 spectrum from a quiescent cloud 



region (f) and the average profile show a decline in in- 
tensity only between —30 and —24 km s~^. The ^^CO 
J=l— >0 spectrum from the peak position of IRAS 16362- 
4845 (e) shows that at least two Gaussian lines - though 
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Fig. 5. 12C0J=2-^1 (top) and i^cQ j^l^O (bottom) emission distribution around IRAS 16362-4845 in the velocity 
range —29.1 to —16.2 km (—30.3 to —17.2 km s~^). The velocity coverage of each map is ~0.5 km s~^ (five velocity 
channels are comprised) and the central velocity is indicated in each panel. Contours go from 3a (1.5 K km s~^) to 
453cr in steps of 45(t for ^^CO J=2^1 and 3a (0.4 K km s^^) to 3G3cr in steps of 30a for "CO J=l->0 . The triangle 
indicates the position of IRAS 16362-4845. 



blended - at -22 and -24 km can be identified. We 
find at this position most intense ^^CO J=2^1 bright- 
ness temperatures of ~30 K. Interestingly, the ^^CO 
J=2^1 spectrum even indicates additional emission in 



the form of a non-gaussian broad wing between -21 and 
-17 km s~^ which is not present in the other spectra. 
Since this emission feature is not well spatially focused 
(see Fig.|S} and only very weak wing emission on the blue 



10 



F. Comeron, N. Schneider, and D. Russeil: Star formation in RCW 108 



side of the spectrum (between -27 and -29 km s^^) is 
found, we assume that it is probably not due to the out- 
flow emission of a YSO. As we will see in the channel 
maps (Fig.O, the -21 and -17 km s~^ emission feature is 
only found at the position of IRAS 16362-4845 and - even 
more prominent - northwest of this source. It may repre- 
sent a second component of the molecular cloud showing a 
close interaction with the HII region: we possibly observe 
swept-up molecular material, arising from clumps eroded 
by ionizing gas from the compact HII region. Since the 
peak of ^^CO J=2-*l and ^^CO J=1^0 emission is not 
found at the position of IRAS 16362-4845, the double-peak 
^^CO J=2^1 emission features can probably not exclu- 
sively be explained by a symmetric expanding shell where 
the compact HII region created a cavity. This scenario is 
suggested by Urquhart et al. H2U(J4|I for explaining their 
CO observations at the position of IRAS source. A more 
complicated geometry with an embedded compact HII re- 
gion and a dense molecular clump at the north-western 
border of the cavity is probably more likely. 

3.2.2. Channel Maps 
Eastern Cloud NGC 6193 

Figure^ shows a series of velocity channels of ^^CO 
J=2^1 (13C0 J=1^0 ) emission around NGC 6193. 
While the emission distribution of ^^CO J=2-^l (top) is 
highly fragmented and dispersed in all velocity ranges, the 
^■^CO J=l— >0 emission (bottom) focusses mainly on a sin- 
gle region close to the pair of O stars in a rather narrow ve- 
locity range (~ — 18 to ~ —15 km s~^). However, there is a 
clear morphological and kinematical coincidence between 
the most prominent structures visible in both molecular 
species. As can be seen in Fig. 4 in Arnal et al. (|2L)L)3I), 
there are molecular clouds further south of our mapped re- 
gion at velocities between —17.5 and —11.2 km s~^ which 
may constitute the remains of an initially more extended 
molecular cloud, now mostly disrupted by the NGC 6193 
cluster. 

The single -24 km s""'^ velocity component visible in 
the spectra of Fig.|3|is prevalent across the whole map and 
looks rather unrelated to the bulk emission of the cloud 
which starts at -20 km s~^ and peaks at —17 km s^^. 
There is a clear velocity shift in ^'^CO J=1^0 emission 
(which is less easily recognized in ^^CO J=2^1 ) in this 
region which becomes even more obvious in a position- 
velocity plot shown in Fig. |H| This cut at constant decli- 
nation (offset 1.'5) shows that the velocity decreases from 
^ —16 km s~^ to ^ —18 km s^^ from east to west. The 
physical interpretation of this gradient is not unambigu- 
ous. If caused by solid-body rotation of the cloud, in which 
the cloud itself is continuous and turns around the O stars 
whose wind/radiation created a hole around them, the ve- 
locity gradient of 1.3 km s^^ across 1' would indicate a pe- 
riod of 1.7x 10^ yr, which is in the same order of magnitude 
like the one for the Rosette Molecular Cloud (3.1x10'' yr. 



BHtz and Thaddeus [TMljl or the Orion Nebula (Kutner 
et al.nWH. 

Expansion powered by the energetic output to the mas- 
sive stars in NGC 6193 provides an alternative explanation 
to the observed structure and kinematics of this cloud, al- 
ready suggested by Phillips et al. (1986) and Yamaguchi 
et al. H1999|l . In this interpretation the arc-shaped peak 
most clearly seen in the —16 and —15.4 km s~^ ^^CO 
J=2^1 channel maps of Fig. 0] would be part of the 
rim of an expanding shell containing HD 150135/150136 
within its contour, whereas the peak close to the posi- 
tions of those stars and having a somewhat more nega- 
tive velocity (~ —17 km s~^) would be part of the shell 
hemisphere located in front of the stars and moving to- 
wards us. Difficulties with this interpretation may be how- 
ever indicated by a comparison of the momentum of the 
shell and the the momentum injected by the winds of the 
O stars on their surroundings during their lifetimes. The 
momentum of the cloud can be estimated by taking the 
velocity difference between the proposed frontal feature 
of the shell and the rim (~ 2 km s"^) and multiply- 
ing it by our mass estimate of the cloud (660 M0; see 
Section l5.2.3|l . For the stellar winds, we adopt typical mass 
loss rates {M ~ 10^^ — 10^^ Mq yr^^), terminal velocities 
{voo 1000 km s~^), and an age of 3.1 x 10^ yr (Vazquez 
& Feinstein 11992(1 . The momentum contained in the ex- 
panding shell turns out to be greater than that injected 
by the stars by about one order of magnitude. The actual 
difference is probably greater, since the currently derived 
660 Mq seem to be only the remnants of a larger, now 
mostly dispersed cloud. On the other hand it is also pos- 
sible to invoke particular density distributions that may 
improve the agreement in our comparison, as would be the 
case if the proposed rim were actually a high density ridge 
expanding at a lower speed than the layer of gas moving 
towards us along the line of sight. Unfortunately, the frag- 
mentary nature of the cloud and the tentative character of 
the interpretation of the structures identified in it prevent 
us from deciding between both outlined possibilities. 

Western Cloud IRAS 16362-4845 

The channel maps of the molecular cloud west of the 
bright rim (Fig. (S)) show two prominent features, firstly 
a sharp gradient in ^^CO J=2^1 and ^^CO J=1^0 line 
intensity at the interface region to the rim, and secondly 
a pronounced peak of both, ^^CO J=2— >1 and ^^CO 
J=1^0 emission close to the location of IRAS 16362-4845 
(marked as a triangle in the figure). Apart of these mor- 
phological elements the western molecular cloud is not as 
strongly fragmented as the eastern one. The main features 
of the molecular gas that we identify near the position 
of the compact HII region are in general coincident with 
those discussed by Urquhart et al. H2()()4|l . However, the 
greater extension of our mapped area and the somewhat 
improved sensitivity in the ^^CO (J = 2 — > 1) line al- 
low us to carry our a more comprehensive discussion on 
the structure and kinematics of the surrounding regions 
as well. 
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Table 1. Excitation conditions of the molecular clouds associated with IRAS 16362-4845 (western cloud) and 
NGC 6193 (eastern cloud). Column 1: Cloud labeling. Cols. 2 and 3: average main beam brightness temperatures 
of the i^CO J=1^0 and ^^CO J=2^1 lines. Col. 4: line center velocity of ^^CO J=1^0 . Col. 5: excitation temper- 
ature derived from ^^CO. Col. 6: optical depth of the ^'^CO line. Col. 7 and 8: ^^CO and H2 column densities. Col. 9: 
total mass. Col. 10: average H2 density assuming a slab with n — N{}i2)/2/r. Col. 11: equivalent radius determined 
by r=(area/7r)°'^ and deconvolved with the beam size. " Taken from Arnal et al. H2(JU3|I . ^ maximum line temperature 
of 12 CO 1^0 Taken from Yamaguchi et al. iflMUji . 
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The HII region NGC 6188 and the molecular cloud 
are in direct contact in the velocity range ~-25 to ~- 
23 km s~^ where we see a sharp gradient in CO inten- 
sity at the interface. Towards higher velocities the inter- 
face region moves westwards and CO emission is only 
found close to the IRAS source. However, low intensity 
emission - best visible in the ^^CO J=2^1 map - re- 
mains across the whole mapping area. ^^CO J=2-^l and 
^■^CO J=l— >0 peak emission is directly associated with 
IRAS 16362-4845 only in the velocity range -26 km s"^ 
to -23.5 km s^^. There, we find ^^CO column densities 
(N(^^CO) = 1.5xlO^'' cm~^) at a moderate optical depth 
(t=0.8). Towards higher velocities, the region of peak 
emission moves north-west and is displaced by ~40-60" 
with regard to the IRAS source. This core region is pro- 
jected over a plateau of lower intensity widespread ^^CO 
J=2— i'l emission. This lower column density material is 
not visible in ^'^CO J=1^0 where the region of peak 
emission is more extended but also more clearly defined. 
We do not see a clear outfiow signature in the ^^CO 
J=2— >1 channel maps (see Sec. l3.2.fl for a discussion of the 
spectra) though a bipolar structure with a blue component 
(-27.4 to -26.3 km s^^) east of IRAS 16362-4845 and a 
more prominent red component (-21.8 to -16.8 km s""'^) on 
the western side can possibly be identified. Higher angu- 
lar resolution should clarify to which extend this emission 
is due to YSO. So far we attribute the redshifted high- 
velocity molecular gas to another molecular cloud frag- 
ment which is eroded by the ionizing gas of the compact 
HII region forming IRAS 16362-4845. 

The morphological relationship between the molecular 
gas mapped in this region and the structures seen in visi- 
ble and near- infrared images is obvious. In addition to the 
noted correspondence between IRAS 16362-4845, its as- 
sociated stellar aggregate revealed by the infrared images 




Fig. 6. Position-velocity cut of ^^CO J=1^0 at constant 
declination (offset 1.'5) close to NGC 6193. Contours go 
from 1 to 7 K km s~^ in steps of 1 K km s~^. 



f Section I3.4|l . and the peak in molecular emission inten- 
sity, the rim nebula corresponds to a steep gradient in 
molecular gas density indicating that the molecular gas is 
dynamically influenced by the presence of the ionization 
front directly ahead of it. The associated shock front prop- 
agating into the cloud is best visible for ^'^CO J=l— *0 at 
w ~ — 24 km s~^. Further west, the CO emission distribu- 
tion is more elongated and looks like streaming away from 
the front at velocities '-^-24.6 to '^-21.3 km s~^. 

3.2.3. Physical properties of the molecular clouds 

Table^gives an overview of the physical conditions of the 
bulk emission of the molecular clouds derived from our 
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i^CO J=2^1 and ^^CO J=1^0 SEST observations and 
in comparison with values from the literature. 

The excitation temperature was calculated by as- 
suming an optically thick ^^CO line so that the ra- 
diative transfer equation simplifies to Tf,x = 11.06 x 
(ln(11.06/(r(i2CO)-)-0.187)-hl))-^ [K] with the line tem- 
perature T(^-^CO) [K] determined with an Gaussian fit to 
the observed line. The ^'^CO line was assumed to be op- 
tically thin and the opacity is then given by t{^^CO) — 
-ln(l - r(i3C0)/(5.289/(exp(5.289/rea;) - 1) - 0.868)) 
with the hue temperature T(^'^CO) [K] equally deter- 
mined with an Gaussian fit. We assume LTE so that the 
excitation temperatures for ^^CO and ^'^CO are equal. 
With the line integrated "CO J^l^O intensity W("CO) 
[K km s~^] over different areas (see Table 1), the ^■^CO 
column density is given by N{^^CO)[cm~'^] = f{T^x) x 
10i5W"i3co (Frerking et al. HW^ with a value of f{Tex) 
of 0.87, 1.06 and 1.87 at an excitation temperature of 10, 
15 and 35 K. The H2 column densities are then calcu- 
lated using the ^^CO column density with N{}i2)[cm~^] — 
4.7 X 105A^(i3CO) [cm-2] (Dickman ITW5|l. The masses 
were determined by M[Mq] = 6.6 x lO-'^^N(}i2)D^A with 
the distance of the cloud D in parsec (1300 pc) and the 
areal extent A in square degrees. The average H2 density is 
evaluated by assuming a slab with a column of the length 
2 X Radius. 

For the western cloud, we distinguish two different ar- 
eas: the cloud core correlated with IRAS 16362-4845 is 
covered by approximately one beam in ^'^CO J=1^0 and 
the whole cloud is defined by the 5a level of the observa- 
tions (essentially all emission visible in Fig.Ol). The latter 
compares best with the literature values. The molecular 
core extends over ~0.24 pc which is approximately the 
same size as for the embedded compact HII region plus 
halo (Sec. 13.1(1 . The extinction given by the H2 column 
density is very high [Ay — 70™), and when compared 
to the extinctions derived from the stellar infrared col- 
ors fSection i3.4|l it suggests that the embedded aggregate 
lies roughly between the edge of the cloud facing us and 
the center of the absorbing column of gas and dust in its 
direction. The volume density is high (9x10* cm^'^) indi- 
cating the existence of a large reservoir of dense molecular 
gas. The cloud core mass is ~200 Mq, which is similar 
to the stellar mass of the whole aggregate fSection 13. 4|) . 
The whole cloud has a mass of ~8G00 Mq which is ap- 
proximately two times larger than the masses derived by 
Arnal et al. ((2()()3|l while the area is two times smaller. 
However, their observations are strongly affected by beam 
dilution (8'. 7 angular resolution) and ^^CO 1^0 was used 
for the mass estimate which may explain the differences in 
masses. Yamaguchi et al. (|1999|l also use -'^^CO J=1^0 as 
a mass tracer and obtain a larger total cloud mass though 
at lower average densities. 

For the eastern cloud, we determined the cloud prop- 
erties for the region of strongest ^^CO J=1^0 emission 
visible in Fig. 0]in the velocity range -18 to -14 km s^^. 
While the optical depth (^0.2) is very similar in compar- 
ison to the western cloud, the ^^CO - and accordingly H2 




Fig. 7. Velocity integrated ^^CO J=2— >1 emission over- 
laid with the limit of the interferometric Ha mapped area. 



- column densities are smaller. The gas is not as dense 
(n ~ 4.5 X 10^ cm-3) as for IRAS 16362-4845 and the 
total mass (660 Mq) is much lower than for the western 
cloud. Also the extinction on the background. Ay — 18™, 
is much lower than in the western cloud. All that, together 
with its more fragmentary morphology, points towards a 
more evolved molecular cloud which constitutes now rem- 
nant material. 

3.3. The kinematics of the ionized gas 

The continuous Ha information over the observed field 
(Figure [7^ allows us to produce maps of the ionized gas 
velocity, line width and intensity (Figure|HJ|. It was always 
possible to fit the observed profiles by a single Gaussian. 

The intensity map (Fig. |S1 bottom) clearly shows two 
major features, the bright rim in the west and the compact 
HII region in the map center. More diffuse emission is 
present in the rest of the field. The compact HII region 
appears displaced from the peak of the CO cloud. However 
we can note that it is located where CO isocontours are 
distorded looking like a cavity from where the ionized gas 
can expand. This is a clear indication of the interaction 
between the HII region and its parental molecular cloud. 

As already mentioned section 3.1., the compact HII re- 
gion is very patchy and slightly elongated in South-East 
direction (Figure IHJ). Integrated over the whole compact 
HII region the Ha profile gives a systemic velocity of 
-21.7 km s^i (FWHM 27.5 km s~i). A detailed analysis 
of the different patches highlights a small velocity varia- 
tion in the same South-East direction with velocities from 
-17.8 to -23.3 km s~^. The brightest clump has a velocity 
of -19 km s"^ (FWHM 29 km s^^). In the direction of the 
opaque dust lane the Ha velocity is also -21.7 km s^^. 
The bulge of CO emission emission at -24.5 km s^^ can 
be associated to the extended emission of the HII region 
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Fig. 9. Full resolution monochromatic Ha image of the 
compact HII region. For some strategic positions, the Vlsr, 
with the line widths in parentheses, are indicated (in 
km s-i). The image is centered on a(2000) = 16''40"0!3, 
(5(2000) -48°52'06" and covers a field of 1'20 x 0'97, 
with North at the top and East to the left. 




Fig. 8. Central velocity (in km s^^), full width at half 
maximum (in km s^^) and intensity distribution (in 
ergs/s/cm^^/sr/km s^^). The overlaid isocontours corre- 
spond to the ^^CO J=2— >1 emission. 



while the Ha patches can be counterparts of the bumps 
around -20 km s~^ seen in the CO profile. This can be 
interpreted as dense clumps externally ionized and pho- 
toevaporated. It is consistent with the conclusion inferred 
from the infrared images where the core of the compact 
HH region is deeply embedded in the molecular cloud and 
only the external, less obscured parts are revealed in the 
visible. 



The line width of the ionized gas (Fig.|Hl middle) shows 
a small and marginally significant increase at the location 
of the sharp edge of the molecular cloud (25 km s~^) while 
towards the West and East the width is around 23 km s~^ . 
This trend is probably due to the presence of the ioniza- 
tion front. Elsewhere in the mapped area and farther from 
the ioniation front, the mean width is 27 km s~^. The line 
width is the combination of two factors, thermal gas mo- 
tions and turbulence. A typical HII region temperature 
of 10000 K gives a thermal FWHM of 21.4 km s'^ If 
we interpret the width broadening by the turbulence this 
leads to a velocity dispersion between 8.5 and 16.5 km s~^. 
Systematic velocity dispersion trends and variations were 
already reported in other HII regions (e.g. Godbout 1 1)91] 
and references therein). 

Finally, a clear large-scale velocity gradient (Fig. |H1 
top) is visible. From velocities of -25 km s~^ in the North- 
East part of the map, we go up to about -11 km s^^ in the 
South-West part with a mean velocity -16.2±4.3 km s^^. 
This gradient points towards NGC 6193 confirming that it 
is the ionizing source. A streaming motion of gas toward 
us can be excluded because such motion would imply a 
velocity gradient with decreasing velocities. The velocities 
at the ionization front up to the compact HII region are in 
good agreement with the CO velocities. The Ha emission 
at these velocities certainly comes from the direct interac- 
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tion between the ionizing flux and the western molecular 
cloud. 

The higher velocities correspond to more diffuse 
and fainter Ha emission except the structure around 
a(2000) = 16''40"08''4, (5(2000) = -48°54'18" (-14 km 
s~^) which is quite intense. The Western cloud does not 
show similar velocities but such velocities are noted for 
the Eastern cloud (~ -15 km s^^). In parallel, Arnal et 
al. H2()()8|l detected molecular clouds further south at ve- 
locities between —17.5 and —11.2 km s^^ which are sup- 
posed to be remains of an initially more extended molec- 
ular cloud, now mostly disrupted by the NGC 6193 clus- 
ter. Hence the Ha emission at velocities between -16 and 
-11 km s^^ is probably the ionized counterpart of this 
initial parental cloud. However due to the limited exten- 
sion of the Ha data it is not possible to arrive to definite 
conclusions. 

3.4. The ionizing aggregate of IRAS 16362-4845 

The stellar component of IRAS 16362-4845 is revealed 
by the JHKs images, where a tight aggregate of stars 
emerges at the center of the compact HII region. As noted 
m Section KT^ the aggregate approximately coincides 
in position with the peak in column density of molecu- 
lar gas, which causes an extinction on the background of 
Ay = 70™ (or Ak = 8™), thus ruhng out that it may be 
actually composed of background stars seen through the 
cloud. 

The color-magnitude and color-color diagrams are use- 
ful tools in examining the membership, rough spectral 
types, and possible existence of infrared excess of the stars 
that appear projected on the nebula. The color-color dia- 
gram of the aggregate at the core of the HII region and its 
surroundings has been discussed by Urquhart et al. (J2004| , 
who have identified several infrared-excess objects based 
on 2MASS data. Our SOFI data reach over two magni- 
tudes deeper, have a better linear resolution by a factor of 
at least two, and oversample the PSF area. Therefore, they 
allow us to obtain a much more detailed view of the stellar 
contents of the aggregate. Both the color-magnitude and 
color-color diagrams are presented in Figures ^1 and ITTI 
for the entire catalog of objects detected in our near in- 
frared mosaic of the region, with the objects projected 
on the central 40" x 40" (0.3 pc x 0.3 pc) of the in- 
frared nebula marked with full circles. The position of the 
members of the aggregate in the color-magnitude diagram 
shows that the brightest stars seem to have late O types, 
and are accompanied by several early B-type stars. The 
tightest concentration occurs in a small cluster precisely 
at the position of IRAS 16362-4845, and occupies a pro- 
jected area approximately 8" across (0.05 pc). Judging 
from their {H — Ks) colors they are reddened by widely 
varying amounts ranging from Ay ~ 8 to Ay — 29 or per- 
haps Ay — 35 (depending on whether the reddest star is 
a true member of the cluster or rather a background star; 
see discussion below), using the Rieke & Lebofsky ((1985)1 




Fig. 10. Color-magnitude diagram of all the objects de- 
tected in the infrared mosaic images. The filled circles in- 
dicate the stars lying projected on the central 40" x 40" 
of the nebula (see also Table ^ . The dashed lines mark 
the expected positions of stars of different spectral types 
(whose unreddened colors and magnitudes are given by the 
nearly-vertical solid line at the distance of Ara OBI) red- 
dened by varying amounts. The range of (H — Kg) colors 
suggest a wide range of extinctions of the stars embedded 
at different depths in the nebula, ranging from Ay ~ 8 to 
Ay ~ 35. 

extinction law and assuming essentially zero {H — Ks) in- 
trinsic colors. The compactness of the aggregate and its 
stellar contents are reminiscent of a Trapezium-like clus- 
ter. Other stars having [H — Ks) colors in that same range 
are visible at some distance from the central cluster but 
still within the confines of the compact HII region. 

3.4.1. Individual members 

The census of stars in the central 0.3 x 0.3 pc^ of the HII 
region is given in Table [3 and plotted in Figure [T^ Some 
of these stars may be either foreground or background, 
and thus unrelated to RCW 108. In particular we con- 
sider Star 13, with blue JHKs colors, as a very likely 
foreground source. In principle we also would be inclined 
to consider as a foreground source Star 8, the one that 
appears brightly in visible images nearly coincident with 
the position of IRAS 16362-4845. However, its {H - Ks) 
is much redder than expected given the (J — H) color, 
hinting at the existence of a near infrared excess indica- 
tive of youth and thus possible aggregate membership. Its 
position in the {H — A'5, Ks) diagram indicates that it is 
unlikely to have a spectral type earlier than AO, as con- 
firmed by its visible spectrum, dominated by Balmer lines 
and by the Call H and K lines. Still near the cluster core, 
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^: Extinction estimate using Rieke & Lebofsky II1985II extinction law and the (H — Ks) color index. The value is given only 
for stars having JHKs measurements whose positions in the color-color diagram indicates no infrared excess. Intrinsic colors 
[H — Ks)a ~ —0.04 and (H — Ks)o = —0.01 are assumed for O stars and stars in the B5-A0 range, respectively. 



Star 5 has H and Ks magnitudes consistent with a late 
O type. We note however that its {H — Ks) color implies 
an extinction Ay = 35 in its direction if entirely due to 
foreground reddening, which would be one of the high- 
est found among the objects in the whole field, although 
only half that inferred from molecular-line observations. 
Unfortunately, the star is too obscured for its J magnitude 
to be measurable and no infrared spectrum is available for 
us to decide on the actual membership of this star. 

Star 12 seems to play the main role in ionizing the 
IRAS 16362-4845 nebula, judging from its central position 
at the core of the nebula and its magnitude. Its photom- 
etry is in agreement with the late-0 spectral type that 
is expected from the excitation characteristics of the neb- 
ula (Section ESI- The fact that IRAS 16362-4845 con- 
tains a tight cluster rather than a single star explains the 
discrepancy between the Lyman continuum flux and the 
luminosity derived from far-infrared data pointed out by 
Straw et al. I)1H87|I . Star 12 is detected in visible-red im- 
ages of the region as a very red object next to Star 8, 
which dominates at visible wavelengths. The visible spec- 
trum of Star 12 shows only a steep rise towards the red but 
no features, and is severely contaminated by the nebula 
at the position of Ha and other strong nebular lines. The 
1.5-2.4 /im spectrum, shown in Figure El is also mostly 
featureless as expected from a late O star. The hint of 



Br7 weakly in absorption must be taken with caution, as 
it may be an artifact due to small-scale structure in bright- 
ness of the nebula affecting the subtraction of the nebular 
lines from the spectrum of the star. However we note that 
the subtraction process removes well the Hel (2.058 /im) 
line, whose intensity in the nebula is similar to that of 
Br7 fSection I3.5|l . leading us to favor the interpretation 
of the Br7 absorption as a real feature in the photosphere 
of the star. Moreover, its strength is consistent with that 
expected for a late-0 spectral type (Hanson et al. 19!l^. 

It is interesting to note the possible existence of two 
faint emission lines in the spectrum of Star 12 at A = 
2.314 fim and A — 2.393 //m. Despite the increased sky 
thermal background on the long-wavelength end of the K 
band the quality of the spectrum of Star 12 seems to be 
sufficiently good at the position of these two lines, espe- 
cially at 2.314 /im. We have attempted several slightly 
differing reductions of the spectra, by choosing different 
background subtraction apertures and by excluding each 
one of the four spectra obtained at a time, in order to 
ascertain the reality of the features, and in all cases we 
have recovered them. The position of the 2.314 /im fea- 
ture matches well that of a feature seen in evolved massive 
stars, such as WRA 751 (Morris et al. 11996(1 . where it is 
nevertheless accompanied by other, much stronger emis- 
sion lines. On the other hand, the good subtraction of all 
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Fig. 11. Color-color diagram of all the objects detected in 
the infrared mosaic images. To avoid degrading the quality 
of the diagram with faint objects whose colors are uncer- 
tain we have limited the sample plotted here to objects 
brighter than Ks = 14.5. The filled circles indicate the 
stars lying projected on the central 40" x 40" of the neb- 
ula (see also Table 
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Fig. 12. ifs-band image of the IRAS 16362-4845 aggre- 
gate, indicating the numbers of the objects whose posi- 
tions and magnitudes are given in Table El 
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Fig. 13. Infrared spectrum of Star 12 (TableEJ, the most 
likely responsible for the ionization of the IRAS 16362- 
4845 nebula. As discussed in the text, the Br7 absorption 
(2.166 /im) seems to be real. This may also be the case for 
the two faint emission lines at 2.314 and 2.393 /im. The 
interval between 1.8 /im and 2.02 fim has been removed 
due to the strong telluric absorption. 
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Fig. 14. Infrared spectrum of Star 16, showing both its 
prominent Br7 emission and the CO bands longwards of 
2.29 /im. The interval between 1.8 /im and 2.02 /im has 
been removed due to the strong telluric absorption. 



other nebular lines suggests an origin in the photosphere 
or the wind of Star 12. Therefore it is unfortunate that we 
cannot offer an explanation for these features, and men- 
tion them here as a possible puzzle should better quality 
observations provide a definitive confirmation of their ex- 
istence. 

Star 16 is almost 2 mag fainter than nearby Star 12 but 
has similar colors, suggesting that both are embedded at 
similar depths in the cloud. Our infrared spectra had the 
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slit oriented so as to contain both stars, and the resulting 
spectrum of Star 16 is shown in Figure Examination 
of the frames containing the spectral trace convincingly 
shows that the prominent Br7 emission in the spectrum 
of Star 16 is related to it, and is not due to a poor sub- 
traction of the nebular lines of the surrounding HII re- 
gion. The spectrum also clearly shows the CO absorption 
bands starting at 2.29 /xm typical of cool atmospheres (e.g. 
Kleinmann & Hall '1986) and suggest that Star 16 might 
actually be a background cool star. However, the simulta- 
neous appearance of both Br7 emission and CO absorp- 
tion is frequent among young stellar objects surrounded by 
significant amounts of circumstellar gas and dust (Greene 
& Lada ll99'^ . which produce a much closer match to the 
spectrum of Star 16. Our data also hint at color variability 
of this source, as may be seen by comparing the spectra 
in Figures E| and with the photometry in Table El 
The imaging observations yield a redder {H — Ks) color 
for Star 16 than for Star 12, while the slopes of their HK- 
band spectra, obtained three years later, show that Star 16 
was slightly bluer than Star 12 in 2003. 

The photometry of Stars 18 and 20 also suggests a 
spectral type earlier than B. Star 18 may indeed be a O- 
type star, but the strong infrared excess derived from its 
JHKs colors, already noted by Straw et al. I|1987|l . in- 
dicates that a substantial fraction of the flux at Ks may 
come from a disk, and that the central star may actually 
be of later type. In either case, the infrared excess can be 
taken as an indicator of its youth and thus of membership 
in the aggregate. The situation of Star 20 is less clear, as 
it displays a slight infrared excess at most, probably in- 
sufficient to explain its brightness at Ks as dominated by 
circumstellar emission. Although it seems to be the intrin- 
sically brightest star of the IRAS 16362-4845 aggregate, 
no obvious effects on the nebula are seen in its immedi- 
ate surroundings, thus leaving open the possibility that 
it could be a background giant instead. We note however 
that the numerous giants seen in the whole field imaged 
in the infrared trace a well-defined reddening vector, from 
which Star 20 is detached by nearly 0.3 mag. The sense 
of the departure from the locus traced by giants is the 
same as observed in large-amplitude variables (e.g. Glass 
et al. ll995|l . This possibility was considered also by Straw 
et al. (|1987(), whose low- resolution spectrum of Star 20 
(= IRS 19 in their list) does show indeed CO absorp- 
tion at 2.3 /im. While considering it unlikely on absolute 
magnitude grounds that Star 20 could be a background 
giant or supergiant, they favored instead actual member- 
ship in the aggregate and considered possibilities such a 
T Tauri star or even a FU Ori object. Further infrared 
spectroscopy of this star should clarify its nature. 

3.4.2. Stellar mass and density of the aggregate 

It is possible to obtain a rough estimate of the stellar 
mass of the IRAS 16362-4845 aggregate by assuming that 
our images are sufficiently deep to record all the members 



with masses above a certain threshold, and then extrap- 
olating the mass function to lower masses to account for 
the fainter, undetected members. Since the distribution 
of colors and magnitudes of stars in the aggregate plot- 
ted in Figure EH suggests that there are practically no 
stars earlier than AO that may have been missed by our 
survey (other than those that may be members of unre- 
solved binary pairs), we will use the mass corresponding 
to that spectral type to derive the scaling factor of the 
mass function, for which we assume a log-normal Miller- 
Scalo form (Miller & Scalo ll979|l . Table El indicates that 
there appear to be approximately 19 stars with spectral 
types AO or earlier, among which we include Star 20 as a 
member but exclude Star 5 as a possible background star. 
The only star discarded as foreground. Star 13, should 
be later than AO if it is actually an unobscured mem- 
ber of the aggregate and thus does not enter our census. 
Similarly, the lightly reddened Star 8 is excluded from the 
counting also due to its confirmed later spectral type. We 
adopt M = 2.9 Mq for a AO main sequence star from 
the compilation by Drilling & Landolt H20()()|l . noting that 
the spectral type-mass relationship may not be extended 
to later spectral types given that stars with lower masses 
take at least a few million years (5 Myr for a 2.5 Mq star; 
D'Antona & Mazitelli lll994jl before reaching the main se- 
quence. 

The mass of the aggregate that we estimate in this 
way is ^ 210 Mq, with a considerable uncertainty due to 
a number of reasons. Our assessment of which stars are 
members and which ones are foreground may be incorrect 
in some cases, thus altering the true census of aggregate 
members. Moreover, as we noted above, we have not made 
any corrections for unresolved binarity, which is known 
to affect a large fraction of massive stars (e.g. Garmany 
et al. I1980|l . Also, the scaling factor of the initial mass 
function is derived from its upper end, which is affected 
by small-number statistics, and some stars may appear 
above the AO limiting line in Figure ^1 because of the ex- 
istence of Ks-ha,nd excess. The evolutionary status of the 
aggregate may also affect our mass estimate in a manner 
similar to that described by Herbig & Terndrup H1986|l 
for the Trapezium cluster in Orion: the high brightness of 
some of the stars that wc detect may be due to their pre- 
main sequence status rather than to their mass, resulting 
in luminosities higher than those of main-sequence stars of 
the same masses, leading to an overestimate of the num- 
ber of massive members. This effect should not affect the 
most massive, 0-type stars like Star 12, whose high tem- 
perature is confirmed by the spectrum of the HII region, 
as such stars have reached already the main sequence by 
the time that they become visible (Palla & Stabler ri990l 
Beech & Mitalas 1994|, but may result in a systematic 
overestimate of the mass of other cluster members and 
therefore of the cluster as a whole. Finally, we note that 
the aggregate seems to contain numerous stars of spectral 
type AO and earlier with {H — Ks) < 1.8, but very few are 
found below this line, and none with H — Ks bluer than 
{H — Ks) ~ 1.2 other than the likely foreground Star 13. 
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No bias in our observations can explain the lack of objects 
in the aggregate in this region of the color-magnitude di- 
agram (which is well crowded for other parts of the field) , 
and we can only explain it as a real absence of moderately 
reddened stars with luminosities below that corresponding 
to a main sequence A star. The apparent absence of stars 
with masses below 2.5 Mq was already noted by Straw 
et al. (^87), who suggested bimodality in the mass func- 
tion as a possible explanation. An alternative intriguing 
possibility to explain these observations might be that at 
the early age of the aggregate only the most massive stars, 
which complete their evolutionary tracks towards the main 
sequence much faster than intermediate-mass and solar- 
type stars, may be sufficiently evolved and emerged from 
their circumstellar envelopes to populate the region of the 
color-magnitude diagram corresponding to moderate ex- 
tinctions. Our observations do not allow us to test this 
hypothesis, but future observations at longer wavelengths 
and high spatial resolution may be able to do it. Both the 
possible overpopulation of the region above the AO main- 
sequence line in Figure 1101 due to the early evolutionary 
stage of the cluster, and the apparent lack of low mass 
stars, may indicate an actual mass of the cluster below the 
210 M0 derived above. In the extreme case that no stars 
less massive than 2.9 M© existed at all in the cluster, and 
assuming that the likely members listed in Table El have 
luminosities near the main-sequence ones, the mass of the 
observed cluster population would amount to ~ 110 Mq. 

It is interesting to compare the mass, contents, and 
extent of the IRAS 16362-4845 aggregate with that of the 
best studied young massive cluster, the Trapezium. The 
mass of the Trapezium can be estimated at ~ 130 Mq 
from the stellar mass density and the approximate ra- 
dius given by Hcrbig & Terndrup 119861 which are respec- 
tively 3000 Mq pc"3 and 0.22 pc. Higher stellar densi- 
ties are obtained when considering only the central re- 
gion of the cluster (McCaughrean & Staufler I1994|l . The 
Trapezium census of the most massive members includes 3 
0-type members among the components of 0^ and 0^ Ori 
(Warren & Hesser I1977|) , some of which are binaries hav- 
ing also high-mass companions (Weigelt et al. 119991 Petr 
et al. lTMHji . The mass that we obtain for the IRAS 16362- 
4845 cluster appears to be somewhat higher, within the 
caveats described above, and the number of O stars is sim- 
ilar, again within the small-number statistics. It should be 
pointed out that the Trapezium does contain a substan- 
tial population of less massive members (McCaughrean 
& Stauffer '1994) such as solar- type stars, low-mass stars, 
and substellar objects, in apparent contrast to the hints 
of a significant lack of stars later than AO that we have 
stressed above. The luminosity function of the Trapezium 
cluster does not suggest any marked deficiencies over any 
mass interval (Muench et al. I2()()2|l . although noticeable 
effects of the most massive stars on the mass function of 
the cluster have been recently suggested by Robberto et 
al. I2004|l . On the other hand, independently of the ac- 
tual stellar content the IRAS 16362-4845 cluster is more 
compact. We estimate a radius of 0.11 pc, only half of 



the Trapezium cluster. If the extrapolation of the mass 
function to lower masses described earlier is approxi- 
mately correct, the mass density is consequently higher 
than in the Trapezium by more than one order of mag- 
nitude, nearly 4 x 10"* Mq pc^"^, which translates into 
riHj = 6.4 X 10^ cm~^ assuming that the mass forming the 
cluster was initially in the form of a single molecular core. 
We note finally that the Trapezium is probably nothing 
more than the nucleus of a much more extended aggre- 
gate containing over 2000 Mq, the Orion Nebula Cluster 
(Hillenbrand & Hartmann I 1 998|l . for which no counterpart 
exists in RCW 108 (See Section EXTJ). 

3.5. Spectroscopy of the HII region 

The long-slit spectroscopy of selected stars in the visible 
and the infrared also provides a cross section spectrum of 
the HII region containing useful information on the phys- 
ical conditions of its different components and on the ion- 
izing stars, complementing that obtained with the obser- 
vations of the stellar aggregate discussed in the previous 
Section. 

The slit position in the visible spectrum contains stars 
8 and 12, and the slit length reaches well up to the edge 
of the western molecular cloud, including a portion of the 
bright rim nebulosity appearing at the erosion interface. 
It thus samples well the compact component of the HII 
region, as well as the rim nebula and the tenuous ionized 
gas foreground to the molecular cloud that pervades the 
whole region. Moreover, when inspecting the spectroscopic 
frames we have appreciated the existence of a concentrated 
knot of emission in the visible appearing where the slit 
position runs closest to Star 17. This knot, whose outskirts 
are included in the slit, can be seen in the y-band image 
(Figure EJ and more clearly in the i?-band image. The 
compact, almost point-like feature in the visible images is 
not coincident with Star 17, which is seen in the infrared 
images only and is actually located 2"0 to the northeast; 
both can be simultaneously seen only in the J-band image 
of Figure 

Our estimates of the density of the ionized gas are 
made on the basis of [SII](A6716)/(A6731) line ratio 
fOsterbrock 1198911 ■ the estimates of the temperature 
of the ionizing radiation arc based on the single-star pho- 
toionization models of Stasinska & Schaerer l|1997|l for a 
gas of solar metallicity and for the density estimated from 
the [SII] line ratios. The translation of an ionizing radia- 
tion temperature into a spectral type of the ionizing star 
is made using the ZAMS models of Schaerer & de Koter 
^1997") . We have estimated the extinctions using either the 
intrinsic H/J/Ha line ratio (foreground nebulosity and rim 
nebula) or the intrinsic Paschen 6/Hq; line ratio (com- 
pact nebula and emission knot near Star 17), taken from 
Osterbrock H1989|l for case B recombination. The extinc- 
tion curve used is that of Cardelli et al. H1989|l . We have 
also considered the empirical extinction curve derived by 
Bautista et al. H1995|l for the Orion nebula, under the ex- 
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Fig. 15. Spectrum of the HII region at the interface where 
the RCW 108 molecular clouds is being eroded by the 
radiation of NGC 6193 and its brightest stars, HD 150135 
and HD 150136. The lower panel uses an expanded vertical 
scale to enhance the visibility of the fainter emission lines. 



pectation that it might better represent the extinction in 
the direction of the compact nebulosity of IRAS 16362- 
4845 as well. However, in the cases where the extinction 
can be simultaneously measured from the Balmer and the 
Paschen line ratios we find highly discrepant values be- 
tween both when using Bautista et al.'s H1995fl curve, 
and a much better agreement when using Cardelli et 
al.'s 1)1989(1 curve, hence our preference for the latter. It 
must be recalled anyhow that the observed spectrum at 
any position is the integration along the line of sight of 
the spectra produced over a range of depths, extinctions, 
and conditions, and that the concept of typical extinction 
towards any ionized component of the region is thus nec- 
essarily an ill-defined one. 

3.5.1. The cloud interface 

The spectrum of the cloud interface, presented in 
Figure [T31 corresponds to that of a lightly obscured neb- 
ula {Ay — 2.2 from the 11(3 /Ha ratio) ionized by mid- 
to-late type stars. The [SII] line ratio indicates a density 
somewhat below 100 cm^'^. The extinction-corrected ra- 
tios of the [OIII](A5007), [SII] (A6716), [Arlll] (A7138), 
Hel (A5876), and [SIII] (A9531) lines with respect to Ha 
are all compatible with an ionizing radiation temperature 
around 40,000 K. The strength of the [Arlll] and Hel fea- 
tures are not consistent with temperatures significantly 
below that value, while [SIII] does not allow for a much 
higher temperature and is in fact better fitted by a tem- 
perature below 40,000 K (see Section r3. 5. 2|l . 

The spectral type corresponding to a temperature of 
40,000 K is 07-08, which must be compared to the spec- 
tral types directly determined for HD 150135/150136, the 
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Fig. 16. Spectrum of the ionized gas foreground to the 
RCW 108 dark cloud. Like in Figure [T51 the lower panel 
uses an expanded vertical scale to enhance the visibility 
of the fainter emission lines. 



two 0-type stars of NGC 6193 that cause the ionization 
of the western cloud. This is in good agreement with the 
spectral type 06/07 given in the literature for HD 150135. 
The spectral type of HD 150136 is often given as 05 in 
the literature (e.g. Hiltner et al. ll969|l . which seems to be 
earlier than suggested by the spectrum of the rim neb- 
ula. The 07 classification given by Whiteoak |1963| for 
this star is in better agreement with the spectrum of the 
nebula. 



3.5.2. The foreground diffuse emission 

The spectrum obtained through the part of the slit that 
runs across the dark areas of the RCW 108 cloud still con- 
tains emission lines similar to a low-intensity spectrum of 
the rim nebula, indicating that the molecular cloud is seen 
through a layer of ionized gas. The [SII](A6716)/(A6731) 
ratio is lower than for the rim nebula, indicating a lower 
density that we estimate at 10 cm^'^ or perhaps even less. 
The intensities of the nebular lines are nearly constant 
across the slit with the exception of a dark patch near the 
rim nebula where they have a noticeably decrease, prob- 
ably due to the existence of a cloud that is embedded in 
this foreground layer and absorbs much of the radiation 
from behind. We have used this local darkening to obtain 
the sky spectrum to be subtracted from the spectrum of 
the foreground emission. This sky-subtracted spectrum is 
shown in Figure [TBI 

The spectrum of the foreground ionized gas is very sim- 
ilar to that of the rim nebula, showing that the source of 
its ionization is probably the HD 150135/150136 pair too. 
The extinction derived from the H/3/Ha ratio. Ay — 2.3, 
is also very similar to the extinction toward the rim neb- 
ula. However we notice that the [SIII] lines in the far red 
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IRAS 16362-4845 (compact HII region) 
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Fig. 17. Spectrum of the compact HII region. Like in 
Figure El the lower panel uses an expanded vertical scale 
to enhance the visibility of the fainter emission lines. 



IRAS 16362-4845 (compact IIII region) 
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Fig. 18. Infrared spectrum of the compact HII region. The 
interval between 1.8 /im and 2.02 /im has been removed 
due to the strong telluric absorption. 

are weaker here. The lower temperature of the ionizing 
radiation that is inferred, around 30,000 K, is in turn in- 
compatible with the line ratios of all the other emission 
lines. Such a discrepancy, although less dramatic, was also 
mentioned concerning the [SIII] lines of the rim nebula, 
which also favored a temperature clearly below 40,000 K. 



3.5.3. The compact HII region 

The extracted spectrum of the compact nebula, shown in 
Figure lTTI corresponds to a spot lying approximately mid- 
way between the closest positions along the slit to stars 17 
and 21 (see Figure [T^ . and is representative of the whole 
compact component. For the background subtraction we 



have chosen a region further to the East that contains the 
typical spectrum of the foreground emission discussed in 
Section 13.5.21 which as pointed out in that Section keeps 
a fairly constant intensity across the area of the dark neb- 
ula. Thus, we expect the spectrum presented in Figure [T7I 
to be essentially uncontaminated by the foreground emis- 



The most obvious differences between the visible spec- 
trum of the compact HII region and the ones of the rim 
nebula and the foreground emission discussed earlier are 
due to extinction, for which we obtain now a value be- 
tween Av = 5.9 (from Paschen 6/Ha) and Ay = 7.8 (from 
H/3/Hq;). The near infrared spectrum, shown in FigureEJ 
is not extracted at exactly the same position as the visible 
one, since the slit was oriented so as to simultaneously con- 
tain stars 12 and 16, and corresponds to a spot about 5" to 
the Northwest of Star 12. It displays the features expected 
in a HII region, with Brackett series lines clearly visible up 
to Br 10, as well as a strong line of Hel at 2.058/zm. The 
Hel(2.058) / Br7 ratio that we measure, 0.63, is near the 
lower edge of the range measured in diverse HII regions. 
This might be a consequence of the softness of the ionizing 
radiation, although Lumsden et al. H2()()3|l caution against 
the use of this ratio to derive the properties of the ionizing 
stars, especially in compact HII regions. 



Extinction, together with the higher density of the 
emitting gas, are the main responsible of the dominance 
of the [SIII] lines on the far red part of the visible spec- 
trum with respect to Ha. The [SII] line ratio is markedly 
different now, and we derive from it an electron density of 
1,500 cm^'^. This is much lower than the 13,800 cm^'^ 
independently derived from radio continuum data by 
Urquhart et al. H2004|l . which may be due to the assump- 
tions implicit in the derivation of electron densities from 
radio continuum data, but also to the possible existence 
of a denser, more deeply embedded component of the HII 
region that does not dominate the visible emission spec- 
trum. Despite the high extinction the H/3 and [OIII] lines 
on the blue part of the spectrum are still visible. H/3 is 
stronger than [OIII] now, arguing for a lower tempera- 
ture of the ionizing radiation than in the cases of the rim 
nebula and the foreground emission. The lower tempera- 
ture is also favored by the other line ratios, particularly 
Hel/Ha which has a good sensitivity to temperature and 
only a mild dependency on density in this range. The tem- 
perature of the ionizing radiation best fitting all the line 
ratios is 35,000 K, and indicates a spectral type later than 
those of HD 150135/150136. Indeed, the temperature is 
in good agreement with a 09 spectral type, which in turn 
matches well the position of Star 12 in the color-magnitude 
diagram discussed in Section The agreement with the 
spectral type estimated by Urquhart et al. H2004|l based 
on the ionizing photon flux is also excellent. 
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Knot near Star #17 
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Fig. 19. Spectrum of the dense knot near Star 17. Like in 
Figure El the lower panel uses an expanded vertical scale 
to enhance the visibility of the fainter emission lines. 



3.5.4. The emission knot near Star 17 

The emission characteristics of the compact knot near 
Star 17 are similar to those of the compact HIT region, 
with the main differences of a stronger extinction {Ay = 
7.7 from the Paschen 6/Ha ratio; the extinction below 
~ 6000 A is too strong for H/3 to be distinguishable). The 
A6731 component of the [SII] pair is now much stronger 
than the A6716 component and the ratio between the two 
is close to the high density saturation value, indicating a 
density in excess of IC* cm~^. The spectrum longwards of 
Ha is rich, making the entire Paschen decrement between 
Paschen 6 and Paschen 12 visible. Taking into account the 
difference in densities inferred for the compact nebula and 
the knot, the relative intensities of the [SII], [Arlll], and 
[SIII] emission lines are all in agreement with the same 
ionizing radiation temperature as derived for the compact 
nebula, T* ~ 35,000 K corresponding to a 09 spectral 
type. Given the proximity of the knot to Star 12 and the 
consistency with its spectral type as deduced from the 
spectrum of the compact nebula, we consider it likely that 
it is a dense clump externally ionized by this star. 

3.6. Star formation across the RCW 108 molecular 
cloud 

Our near infrared observations, covering a projected area 
of nearly 5 pc x 5 pc, allow us to investigate the traces 
of recent and ongoing star formation well beyond the im- 
mediate vicinity of IRAS 16362-4845, at a scale that is 
intermediate between the star formation directly associ- 
ated to the compact HII region and the large-scale, low 
resolution surveys of Yamaguchi et al. H1999|l and Arnal 
et al. H2003fl that encompass the entire RCW 108 complex 
and most of the Ara OBI association. While IRAS 16362- 



4845 is an obvious site of current massive star formation, 
other sites in the same cloud being obscured by similar 
amounts and forming lower-mass stars may exist as well 
without showing so obvious signposts. One of this sites 
is CD-48 11039, already mentioned in Section ITTl which 
is recognized by the reflection nebula that it illuminates. 
Other young stellar objects of even lower mass may pass 
completely unnoticed. 

To try to identify young lower mass stars in the region, 
we use the position in the (J — H), {H — Ks) diagram 
as a telltale diagnostic for the existence of hot circum- 
stellar disks remnant from their formation, a very com- 
mon signature of youth among intermediate-mass, pre- 
main sequence stars and a frequently used approach to 
the identification of distributed star formation in molecu- 
lar clouds (e.g. Li et al. ll997[ Massi et al. [^OOOt Brandner 
et al.EOnil Jiang et al. 12002 see also Lada & Lada lM75l 
for a review) . To produce a reliable identification of near 
infrared excesses we have considered only stars brighter 
than Ks = 14.5. The luminosities to which this limit cor- 
responds is obviously highly dependent on the amount of 
light reprocessed by the circumstellar material into the Ks 
band and on how deeply embedded in the cloud the star is. 
As a representative number, a Ks = 14.5 star 1 Myr old 
having a foreground obscuration of Ay — 20 {Ak = 2.2) 
as typically found in the IRAS 16362-4845 aggregate, and 
in which half of the emitted flux at Ks comes from repro- 
cessed light would have a luminosity L ~ 1.2 Lq^, which 
corresponds to a mass of 0.8 M© at an age of 1 Myr. 

To assess the amount of infrared excess we use the 
reddening-free quantity 

Q = {J-H)-1.70{H-Ks), 

which measures the separation between the position in 
the color-color diagram of a star with colors (J — H), 
{H — Ks), and a reddening vector that traces the Rieke 
& Lebofsky ('1985') extinction curve having its origin at 
the intrinsic colors of a AOV star. Background red giants, 
which align along a narrow strip running above this red- 
dening vector in the {J — H), {H — Ks) diagram, have 
Q > 0, while early-type stars with no infrared excess 
cluster around Q — 0. Stars with Q < can be either 
late-type M dwarfs (Bessell & Brett I1988|l or stars with 
infrared excess. The former are far too faint to be de- 
tected in our observations at the distance of RCW 108, 
and the areal density of field late-M dwarfs at the limiting 
magnitude Ks = 14.5 is too low for them to appear in sig- 
nificant numbers in the imaged field (Reid et al. |^02). A 
third possibility is that they may be binaries unrelated to 
RCW 108 consisting of pairs with widely different infrared 
colors, which should be a rare occurrence as well. 

Out of the 4,365 stars brighter than Ks = 14.5 for 
which complete JHKs photometry is available, we have 



We take a bolometric correction in the Ks band BCk = 
2.1 as a representative value for a 4,400 K pre-main sequence 
star (Kenyon & Hartmann 119951 in the convective part of its 
evolutionary track (D'Antona & Mazitelli ll994ll 
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Fig. 20. Ks-hand mosaic centered on the position of IRAS 16362-4845, containing most of the RCW 108 cloud and 
parts of its surrounding area. The circles mark the position of stars brighter than Ks = 14.5 displaying infrared excess 
emission according to the reddening- free Q parameter defined in Section ITHl 



found 87 satisfying the Q < —0.10 criterion that we im- 
pose (rather arbitrarily as far as the absolute value is con- 
cerned) as a threshold defining the stars suspected to dis- 
play infrared excess. Their spatial distribution is plotted 
in Figure EOl superimposed on a iiTs-band image of the re- 
gion centered on IRAS 16362-4845 and where the extent of 
the dense parts of the RCW 108 clouds is well delineated 
by the areal density of stars. 

The distribution of stars with infrared excess in 
RCW 108 is clearly non-random. The greatest concentra- 



tion takes place in the densest region of the cloud, near 
IRAS 16362-4845, where the overall density of stars is 
lowest due to the presence of the obscuring cloud. The 
infrared excess stars are not uniformly scattered in the 
cloud, but are instead distributed along a belt that con- 
tains IRAS 16362-4845 and runs from northeast to south- 
west. The northeastern side is particularly interesting, as 
it runs roughly parallel to the rim nebula in the inner 
side of the cloud, suggesting that the stars in this re- 
gion may have their origin in the action of the ioniza- 
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tion front encountering molecular gas, perhaps producing 
the implosion of preexisting cores by the increased ex- 
ternal pressure. We note in particular a clear contrast in 
density of near-infrared excess stars as we move along an 
East- West strip immediately North of IRAS 16362-4845, 
in which the presence of such objects ends about 1' East 
of the HII region without any corresponding decrease in 
molecular gas column density as we proceed further to 
the West. This suggests that, despite the similar physi- 
cal conditions of the molecular gas along this strip, star 
formation has taken place only in the parts closest to 
the ionization front. It is interesting to note that a sim- 
ilar conclusion is independently reached by Urquhart et 
al. (20® based on the location and expected coevality 
of the three thermal sources detected in MSX observa- 
tions (one of which is IRAS 16362-4845), which they clas- 
sify as possible ultracompact HII regions corresponding 
to the earliest stage of massive stellar evolution, at ages 
< 10^ years. One of our infrared-excess stars in this zone, 
at a(2000) = 16''40™11''7, (5(2000) = -48°48'58", is sur- 
rounded by a tiny reflection nebula seen only in the J- 
band image. Some stars with infrared excess are seen also 
outside this area, where the ionization front has already 
destroyed the cloud, and they may be the result of past 
episodes of star formation when the cloud extended fur- 
ther towards NGC 6193. Indeed, assuming that the freshly 
ionized gas can escape freely leaving an essentially clear 
line of sight between the ionizing stars and the rim nebula, 
taking our derived density of nj^^ — 3.9 x 10^ cm^^ for 
the molecular cloud and assuming that the true distance 
r from the ionizing stars to the rim nebula is roughly the 
same as the projected distance, we obtain a propagation 
velocity v of the erosion front 



where S'» = 1.1 x 10^^ s~^ is the flux shortwards of 
the Lyman continuum of a 07 star (Schaerer & de 
Koter [T^97|) . This is an upper limit to the velocity, since 
the low density column of ionized gas between the stars 
and the rim nebula must be kept ionized thus decreasing 
the ionizing flux that reaches the molecular cloud, and the 
true distance between the ionizing stars and the rim of the 
molecular cloud is likely to be somewhat larger than the 
projected distance'^. Nevertheless, the time needed for the 
ionization front to travel at this speed from the position 
of the easternmost infrared excess stars to the present po- 
sition of the rim of the cloud is only ~ 2.7 x 10^ years (of 
the same order of the estimate by Urquhart et al. r2004) 
based on a different set of hypotheses), or over one order of 
magnitude less than the typical disk dispersal timescales 
of solar-type stars (Strom et al. ll993|l . so even a decrease 
of the propagation velocity of the erosion front by the 

* Based on their interpretation of infrared emission produced 
in the photodissociation region at the rim nebula/molecular 
cloud interface, Urquhart et al. ( 2004) suggest that NGC 6193 
lies actually somewhat behind the cloud. 



same factor would still be consistent with the presence of 
stars with infrared excess far beyond the present edge of 
the molecular cloud. Overall, we find 39 stars (45% of the 
total) on the side of the cloud facing NGC 6193 and be- 
yond the cloud in the same direction, and thus possibly 
having a triggered origin. 

It is interesting to note however that numerous infrared 
excess sources are found in the Southwestern quadrant of 
Figure|2n| The Ha image (Figure^]) shows that this region 
is still fairly opaque, despite the higher density of stars 
shown by the infrared images. Also CO maps such as those 
shown in Figure|Slshow streamers of molecular gas extend- 
ing in this direction. Based on their location it seems more 
difficult to link the existence of these stars, and in gen- 
eral all the stars located to the west of IRAS 16362-4845, 
with the action of the HD 150135/150136 pair. No other 
0-type stars have been identified near the Southwestern 
edge of the RCW 108 cloud that may cause effects similar 
to those of HD 150136/150136. The near-infrared stars 
near the Southwestern edge of the cloud may thus rep- 
resent star formation not being triggered by an exter- 
nal cause. Evidence for the coexistence of different modes 
(triggered/non-triggered) of star formation has been also 
reported in other massive star formation regions (Jiang et 
al. 1200 21. However, the fact that such stars are also non- 
randomly distributed across the cloud, and that there is 
no distinct spatial separation between them and the stars 
for which we proposed a triggered origin, puts a note of 
caution in our interpretation of the latter, whose location 
on the side facing NGC 6193 may also be fortuitous. 

3.6.1. Star formation efficiency in the western cloud 

The derivation of the stellar mass of the RCW 108 complex 
on the basis of the members identified via their infrared 
excess would require determining their masses as well as 
the fraction of cloud members that display such infrared 
excess. Unfortunately, the currently available data prevent 
this exercise. However, it is possible to use our estimate of 
the total mass of the cloud (Section I3.2.3|l together with 
indirect arguments to show that the western cloud has in 
any case a low star formation efficiency (defined as SFE — 
M,/(Mgas + A/*), where Af* and Mgas are respectively the 
masses in stellar and gaseous form) as compared to typical 
evolved giant molecular clouds. The stellar mass estimate 
for the IRAS 16362-4849 aggregate, 210 Mq, represents 
only 2.6 % of the mass in molecular form (Section l3.2.3(l . 
well below the ~ 10 % of gas that is typically turned 
into stars in a giant molecular cloud at the end of its life 
(WiUiams & McKee [TM7|l . The fraction is even smaller if 
the apparent deficiency of low mass stars in the aggregate 
hinted at in Section POI is real. In order to raise the star 
forming efhciency to 10 % one thus needs to invoke the 
existence of stars amounting to nearly 700 scattered 
across the cloud outside the boundaries of the aggregate. If 
this component obeyed a Miller-Scalo mass function, the 
expected number of stars more massive than 3 Mq would 
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be around 50. Our observations indicate that the number 
is far smaller than that, with CD-48 11039 being the only 
star outside the aggregate with a mass in this range. The 
low star-to-gas mass ratio, together with the vast amount 
of molecular gas available to star formation and its high 
volume density, thus suggest that most of star formation 
in the ROW 108 still has to take place. It is interesting 
to speculate that, should star formation proceed across 
RCW 108 in the future until reaching a final '--^ 10 % 
efficiency, the result may then be similar to the extended 
Orion Nebula Cluster. The role played by IRAS 16362- 
4849 as the core of that future RCW 108 aggregate would 
then be similar to the role that the Trapezium plays in the 
Orion Nebula Cluster nowadays. In this respect, RCW 108 
and the Orion nebula complex might thus be regarded as 
similar structures at different stages of their evolution. 

4. Summary and conclusions 

In this paper we have presented a collection of visible, 

near infrared and millimeter observations of the RCW 108 
region and its associated rim nebula, focusing on its ac- 
tive massive star forming site, IRAS 16362-4845, and on 
the molecular gas presumably associated to the cluster 
NGC 6193. IRAS 16362-4845 is a compact HII region 
whose actual morphology and stellar contents are best re- 
vealed by near infrared observations. At the core of the HII 
region lies a Trapezium-like compact cluster dominated 
by a late 0-type star, surrounded by a looser aggregate 
probably including between one and three other late O- 
type stars and about 16 B-type stars. The near-infrared 
color-magnitude diagram of the aggregate gives an esti- 
mate of its mass of 210 M0. The cluster seems to lack 
moderately reddened stars later than AO. 

We have discussed the spectrum of the ionized gas in 
the region and the physical conditions of the molecular gas 
as inferred from ^^CO J=2^1 and ^'^CO J=1^0 map- 
ping. The compact HII region is dense (--^ 1500 cm^'^), 
although much less dense than the surrounding molecular 
gas as inferred from the molecular-line observations in the 
same direction. It spectrum suggests a spectral type 09 
for the ionizing star, in good agreement with the position 
of the brightest star of its cluster in the color-magnitude 
diagram. A compact knot located near the main ionizing 
star, whose density is derived to be above 10"* cm ~^ and 
whose lines indicate ionization by a spectrum similar to 
that of the stars ionizing the compact nebula, may thus be 
externally ionized by those same stars. Streaming motions 
are revealed by the interferometric Ha observations both 
in the compact HII region and near the molecular cloud 
interface, which can be interpreted in terms of photoevap- 
oration of the molecular gas. 

We obtain masses of 8000 Mq and 660 Mq for 
the molecular clouds associated to RCW 108 and to 
NGC 6193, respectively. The extinction on the back- 
ground at the position of IRAS 16362-4845 reaches up to 
Ay = 70 mag, while the peak extinction produced by the 
cloud associated to NGC 6193 is much lower, Ay — 18. We 



consider explanations to the kinematics of the NGC 6193 
cloud based on rotation and on expansion powered by the 
stars in the cluster, although we do not find conclusive ev- 
idence for either interpretation. We attribute broad wing 
CO emission (redshifted with regard to the bulk emission 
of the molecular cloud) northwest of IRAS 16362-4845 to 
erosion of clumps due to the ionizing gas of the comapct 
HII region. A more complicated scenario including an ex- 
panding shell and outflow emission from YSO can not be 
excluded but needs further observational investigation. 

Using near infrared excess as a way to identify young 
stars still surrounded by circumstellar dust, we find that 
their distribution is inhomogeneous: not surprisingly, the 
molecular cloud harbors most of them, mainly in the sur- 
roundings of IRAS 16362-4845. Many are distributed in a 
broad band near the edge of the molecular cloud and par- 
allel to it, suggesting that their formation may have been 
triggered by the progress of the ionization front traced 
by the rim nebula into the molecular cloud. Other in- 
frared excess stars appear outside the cloud on the side 
facing NGC 6193, suggesting that they formed in ancient 
regions of the RCW 108 molecular cloud that have been 
eroded away. However, numerous infrared excess sources 
appear also on the side of the RCW 108 cloud opposite 
to NGC 6193, where there is no obvious triggering can- 
didate. Despite the evidence for widespread ongoing star 
formation we find a low overall star formation efficiency 
that, together with the physical conditions and mass of 
the cloud, suggests that most of the star formation in 
RCW 108 still has to take place. Based on this, we spec- 
ulate that RCW 108 may represent a structure similar to 
the Orion nebula complex observed at an earlier stage of 
evolution. 

The observations presented and discussed here, while 
not being an exhaustive survey of the different compo- 
nents of RCW 108 and the gas associated to NGC 6193, 
provide some useful elements towards a complete picture 
of a complex that features an emerged cluster with mas- 
sive stars and the last remnants of its parental cloud, lying 
on the sky side-to-side to a young, embedded massive star 
forming region in earliest stages of interaction between a 
newly formed stellar aggregate and its surrounding gas. 
At a larger scale, we can also investigate possible evidence 
of star formation in a molecular cloud being triggered by 
an external action, namely the ultraviolet radiation of the 
hottest stars of NGC 6193. The identification of intermedi- 
ate mass stars with near infrared excesses across the cloud 
gives us a possibility of carrying out follow-up observations 
aimed at deriving their approximate ages, which might in 
turn provide evidence for a wave of star formation run- 
ning across the RCW 108 cloud. Finally, a detailed study 
of the RCW 108 region can complement observations at 
larger scales comprising the entire Ara OBI association, 
to which the RCW 108 complex belongs, thus yielding 
so far elusive observational views on the interplay of star 
forming processes at widely different length scales. 
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